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It is shown in tliis paper that the transition kernel corresponding 
a, spatially inhomogene' 
lower Gaussian estimates. 



to a spatially inhomogeneous random walk on Z'' admits upper and 



1. Introduction. We shall consider in this paper spatially inhomogeneous 



l> 

(N 

S_ . 

, random walks (5j)jgN with bounded symmetric increments in Z^. More 

1^ ' precisely let T = —T C Z'^ be a symmetric finite subset of Z*^ and let 

^ i r — > [0, 1] such that 

^7r(x,e) = l, 7r(x, e) = 7r(x, — e), e€T,xGZ'^. 

eer 

. Then we let (S'j)jgN be the Markov chain defined by 



P[S'j+i = x + e/ /Sj =x\ =7r(x,e), e£T,x£ Z'^J = 0,1, . . . . 



I To avoid unnecessary complications we shall assume that T contains and all 

■ unit vectors in Z'^, that is, all e with |e| = 1 where | • | denotes the Euclidean 

\^ , norm. Furthermore we shall impose the following ellipticity condition: 



(1.1) 7r(x,e)>Q, x£Z'^,e£r, 

for some q > 0. It must be emphasized that the random walk (5j)jgN is not 
necessarily reversible. 
We shall denote by 

(1.2) pn{x,y) = P,\Sn = y], n = \,1,...,x,y^Z\ 
the transition kernel corresponding to the chain (S'j)jgN and by L the cor- 



\ responding generator, that is, the difference operator defined by 

(1.3) Lfix) = J2 e)(/(x + e) - f{x)), / : Z'^ ^ R. 



eer 
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2 S. MUSTAPHA 

We shall prove that there exists a unique (up to a multiplicative constant) 
positive solution M(-) of the adjoint equation 

(1.4) L*M(x) =^7r(x-e,e)M(x-e) -M(x) = 0, a; G Z"^, 
eer 

globally defined on TJ^ (cf. Section 3.2). We shall denote 
V{x,r)= M{z), xeZ'^,r>0, 

where Br{x) = {y G Z'', \y — x\ < r},x G Z'', r > 0. 

Theorem 1. Let {Sj)j^j<s be as above. Letpn{x,y), x,y £ Z'^, n = 1,2, . . . , 
denote the corresponding transition kernel. Then there exists C > 0, depend- 
ing only on d, T and a, such that 

CMjy) f_\x-y[^' 
P-^'^^y)^ (V{x,V^)V{y,V^)y/^''''^[ Cn 

(1-5) 

n > l,x,y G Z"*, 

Mjy) ( C\x-y\^ 

C{V[x,y/n)V{y,y/n))^/^ \ n 

(1-6) 

n>l,x,y£Z,\x — y\ <n/C. 

The following comments may be helpful in placing the above theorem in 
its proper perspective. 

(i) It will be shown in Section 3 that the volume function V{x, r) satisfies 
the doubling property 

(1.7) V{x,r)<CV{x,2r), xGZ'^,r>0. 

The volume factor (y(x, y^)y(y, v^))^/^ in (1.5), (1.6) can therefore be 
replaced by V{x.,^/n). 

ill) In the reversible case. Theorem 1 is an immediate consequence of 
Delmotte's work (cf. [5]). Delmotte proved the equivalence of the upper and 
lower Gaussian estimates to the volume doubling property (1.7) plus the 
Poincare inequality for reversible Markov chains with bounded increments 
on graphs. His approach relies on a clever adaptation of the Moser iteration 
process. The reversibility property p{x,y)m{x) = p{y,x)m{y) verified by the 
kernel of the chain {Sj)j£]si and its invariant measure m plays a crucial role 
in [5] (cf. also [2]). 
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(iii) Reversible Markov chains on Z*^ are the discrete analogues of diffu- 
sions generated by second-order differential operators in divergence form and 
the inhomogeneous walks can be considered as the analogues of second-order 
operators in nondivergence form (cf. [12], Table 1, page 78). The first two- 
sided Gaussian bound for fundamental solutions of parabolic equations in 
divergence form with measurable coefficients is due to Aronson (cf. [1]). For 
operators in nondivergence form such upper and lower Gaussian estimates 
were proved only recently by Escauriaza in [6]. 

(iv) In Aronson's work the parabolic Harnack inequality is used to obtain 
the Gaussian lower bound (cf. [1]). In fact both upper and lower bound for 
the heat kernel can easily be deduced from the parabolic Harnack principle 
(cf. [19]). Conversely it is shown in [10] that the two-sided Gaussian bound 
implies the Harnack inequality. Saloff-Coste showed in [16] that the parabolic 
Harnack inequality (or the two-sided Gaussian bound) for a divergence-form 
second-order operator (or for the Laplace-Beltrami operator on a Rieman- 
nian manifold) is equivalent to a family of Poincare type inequalities for 
balls and the doubling property (cf. also [11]). The results of [5] are the 
discrete counterpart of [16]. It will be interesting to discuss this type of 
equivalence for both nondivergence differential operators and nonreversible 
random walks. 

A general outline of the paper is as follows. Section 2 collects the main 
potential theoretic properties of spatially inhomogeneous random walks on 
TJ^ . The two new results of this section (Theorems 4 and 5) are of indepen- 
dent interest and are proved in Section 5. In Section 3 we define the concept 
of normalized adjoint solution adapted to spatially inhomogeneous random 
walks and we prove that adjoint solutions verify a parabolic Harnack prin- 
ciple. This Harnack principle is used in Section 4 to deduce the Gaussian 
estimates of Theorem 1. 

2. Potential theory. Let AdTJ^ denote a bounded domain (i.e., a finite 
connected set of vertices in Z'^). We let 

dA = {x E ^ X = z + e, for some z £ A and e G T}, 

r being as in the previous section, and 

A = AudA. 

Let B = Ax {a< k<b} CZ'^ X Z where AcZ'^ and where a<b£Z. We 
let 

diB= U dAx{k}, 

a<k<b 



dpB = diBU{Ax{a}), 
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and B = BU dpB. dpB is the parabolic boundary of B and diB is its lateral 
boundary. We say that u : A — > R is harmonic va. AdTJ^ if 

Lu(x) = ^ 7r(x, e)(u{x + e) — u(x)) = 0, x G A. 
eer 

Let B = Ax{a<k<b}cZ'^xZ and u : S — > R. We say that u is caloric 
in B if 

Cu{x, k) = ^ 7r(x, e)u{x + e^k) — u{x, k + 1) =0, 
eer 

{x,k) £ A X {a<k < b}. 

The following maximum principle is immediate. 

Theorem 2 (Maximum principle). Let B = A x {a < k < b} C Z'^ x Z, 
where a < 6 G Z and A is a bounded domain in Z'^ and let u: B — > R such 
that Cu = in B and u>0 on dpB. Then u>0 in B. 

The following theorem (cf. [13]) is a random walk version of a well-known 
and fundamental result in the potential theory of second-order equations in 
nondivergence form (for an elliptic version cf. [14]). 



Theorem 3 (Parabolic Harnack principle). Let u be a nonnegative caloric 
function in B2r{y) x {s — 4r^ <k<s}, {y,s) e Z'^ x Z, r > 1. Then 

sup{ii(2;, fc); X G Br{y),s — 3r^ < k < s — 2r^} 

(2.1) 

< C mf{u{x,k);x £ Br{y), s — r <k<s}, 
where C = C{d, a, T) > 0. 



In the proof of Theorem 1, together with the previous results we need the 
following estimates which describe the boundary behavior of nonnegative 
caloric functions. Let yQ G Z'^, let iio > and let Q = i?R(,(yo)- 

Let Q = X Z, let Y = {y,s)£dnxZ and let c < r < Ro/2 where c> 
denotes a constant times diam(r). We shall denote 

C,(y) = Briy) x{s-r^<k<s}, Q,(y) = QnCr{Y), 
Yr = {yr,s + 2[r]2), Y, = {yr,s - 2[rf), 

where yr £^ satisfies \yr — {R — r/2) | < 1 and [r] denotes the greatest 
integer < r. 
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Theorem 4 (Boundary Harnack principle). Let Y = {y,s) e dQ x Z. 
Let c<r<Ro/K where K > is large enough. Assume that u and v are 
two nonnegative caloric functions in Q H (B^Kriv) x {s — 9K^r^ < k < s + 
9i^V2}) andu = on {dQ x Z) n {B2Kr{y) x {s - AK'^r'^ <k<s + AK'^r'^}) . 
Then 

(2-2 sup -<C—— — -, 

where C = C{d, a, T) > 0. 

Theorem 5 (Backward Harnack principle). Let u be a nonnegative caloric 
function in Br{yo) xN (yo ^ Z'^ , r > large enough) vanishing on dBj.{yo) x 
N. Then 

(2.3) u{x,k + 2[rf)<Cu{x,k), {x,k) e Br{yo) x {r^ < k <3r'^}, 
where C = C{d,a,T) > 0. 

The proofs of the estimates (2.2) and (2.3) which follow strongly the proofs 
of the corresponding facts about the boundary behavior of nonnegative solu- 
tions of second-order equations in nondivergence form (cf. [3, 4, 7, 8, 9, 15]) 
are given in Section 5. 

We shall use throughout the usual convention / ~ g to indicate that 
<f/g<C for an appropriate constant C > and C, c are used to 
denote different positive constants which depend only on d, a, diam(r). 

3. The adjoint Harnack principle. Let D CZ"^ denote a bounded domain 
and a < 6 G Z. We say that v = v{x, t) : Z) x {a < t < 6} — > R is a parabolic 
adjoint solution of L m D x {a<t<b}, if f satisfies the equation 

v{y, t + l) - v{y, t) = L*v{y, t), t = a, . . . ,b - l,y G D, 

where L* is defined as in (1.4). Let m(-) be a fixed positive adjoint solution 
for L in D [i.e., m : D ^ R, m{x) > 0, Vx G -D, L*m = in Z)]. For instance, 
if Z) C i?ro (0) liss in the Euclidean ball centered at and of radius ro > 
large enough, we can set m{x) = G{x*,x) where x* G B^roiO) \B3roiO), with 
G{-,-) being the Green function of (S'j)jgN in the ball B^roiO). Let u be a 
parabolic adjoint solution in D x {a < t < b}; the function 

v{y,t) = ^^^, {y,t)eDx{a<t<b}, 

is called a normalized parabolic adjoint solution L in D x {a <t <b} (cf. 
[4]). 
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Theorem 6. Suppose that v is a nonnegative normalized adjoint solu- 
tion for L in Br{yo) x N, where yo G Z'^ and r > large enough. Then there 
exists a constant C > depending only on d, a, T such that 

sup{v{y,s);y e B,,f2{yo),r'^ <s< 2r'^} 

(3.1) 

< Cmf{v{y, s);ye 5,,/2(?/o), Sr^ < s < Ar^}. 

Proof. Let Br{yo) = and let qt{-, ■) denote the Green function of C 
in Br- An easy induction on t gives the following representation formula for 
the parabolic normalized adjoint solutions: 

v[y,t)= > m{x)v{x,0) 

(3.2) m{x)v{x,s) Y: 7r(x,e)^i^^^i^f±^, 

yG -Br/2,* =1,2,..., 

where = {e £T /x + e G Br} ■ On the other hand, let us observe that if we 
extend qsi-,y) by qs{-,y) = 0, s < 0, in an appropriate neighborhood of dB^ 
and use the boundary Harnack principle (2.2) to compare ui{x,t) = qt{x,y) 
and U2{x,t) = qt+ior^i^iy)i combined with the backward Harnack principle 

(2.3) , we deduce then that 

qs{x,y) <Cqs+r'2{x,y), dist(x,95r) < cr, 

(3.3) 

y G Br/2,0 < s < 2r2. 

Let now yi, y2 G Br/2 and <ti< Ir"^, 3r^ <t2< Ar^. By (3.2), (3.3) and 
Theorem 5 we have 

v{yuti)< E m{x)v{x,0)^i^^j^ 

(3.4) 

q2ix + e yi) 

+ X! m{x)v{x,s)Y Tiix^e)— — , ' \ 

qt2{x,y2) 



(3.5) 



i^(y2,i2) > TO'(x)?;(x,0) 

,1 / N / .qr^{x + e,y2) 

+ 7?Z^ 1^ m{x)v{x,s) Tr{x,e) 
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A simple use of the boundary Harnack principle [combined with (2.1) and 
(2.3)] allows us to deduce then 

v{yi,ti) < C 7v{y2,t2). 

m{yi) qr2{yo,y2) 

Thus, to prove (3.1) it suffices to show that 

foa\ QAyQ,yi) qr^{yo,y2) ^„ 

The estimate (3.6) is a consequence of (3.4), (3.5) and the fact that 1 is 
a normalized parabolic adjoint solution. Indeed, let 'tp{x) = G(x,yo) be the 
Green function of L in Br with pole at yo- Let A£ Br fixed with | A — yo| ~ 
r/4. By the boundary Harnack principle and the backward Harnack principle 
we have 

QrnX, y) ~ ^^^^ , xeBr\ 53r/4,y e Br/2- 

If we apply now (3.4) and (3.5) to u = 1, we deduce then 
+ Cr'^ m{x) Y ■K{x,e)qr2{yo,yi 



m{yi)<C m{x)——qr2{yo,yi)+C ^ m{x)qr2{yo,yi] 



ip{x + e) 

1 ^ .i^ix) 1 



m{y2)>— J2 'nT'{x)—7j7qr2{yo,y2) + ^ J2 rn{^)Qr^{yo,y2) 

xe-Br--B3r/4 ^'G-B3r/4 

+ 77 m(x) 2^ 7r(x,e)g^2(7/o,7/2) .... , 
and these two inequalities imply (3.6). □ 



3.1. The doubling property for the adjoint solutions. We start with a 
doubling property for the Green functions. 

Proposition 1. Let R be large enough and let xq € Z'^. Let G^{-,-) 
denote the Green function o/(S'j)jgN in the ball B4r{xo). There exists then 
a constant C > ( independent of xq and R ) such that 

Y G^(x,y)<C Y G''ix,y), xeB^R{xo),l<r<R/2. 

y&B2r{xo) y^Br{xo) 

(3.7) 
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This proposition is a consequence of the parabohc Harnack principle and 
the following lemma. 

Lemma 1. Let R be large enough and let xq G Z'^. Let hf {x,y), 
t = 0, 1, . . . , 3;, y G -64^(3:0), denote the heat kernel of {Sj)j£]si in the ball 
B4,r{xq). Then there exists c> (independent of xq and R) such that 

(3.8) inf V hf{z,y)>c, 1 < s < 1 < r < i?/2. 

Indeed, by the parabolic Harnack principle we have 

C inf hg,{z,y)> sup Yl ''M^^v) 

z^BM^o) ^^^^^^^^ zeB,A^o)y^BAx,) 

> inf J2 h%iz,y)>c. 

We have then 

E hf^^^.,{x,y)> J2 E hf{x,u)hg,{u,y) 



> 



E i E hg,{u,y)\hf{x,u) 



U&B2r{xo) 

and this implies the doubling property for the heat kernel with the time 
shifted t^t + 2r2. If we sum on t we deduce that 

00 

Y G^ix,y)<c Y E ^t i^^y) 

y&B2r(xo) y&Br{xo)t=2'P 

<C Y G''{x,y). 

y&Br(xo) 

The adjoint Harnack principle (3.1) and Proposition 1 give 

Theorem 7. Let L*m = 0, m>0, in B^riz), z e Z'^ and r > large 
enough. Then 

(3.9) Y Hy)<c Y ^(y)^ 

y&B2r{z) y&Br{z) 

where C = C[d, a, T) . 
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Indeed let x* G i?6r(-z) \ ^5r(^)- Let G{-,-) denote the Green function of 
(S'j)jgN in the bah Bjr{z). We have 

The second inequahty fohows from the adjoint Harnack principle and the 
third one from the doubling property (3.7). 

Proof of Lemma 1. Let 

U{x,s)= hf{x,y), xeB4r{xo),s>0. 

Let V{x,s) be the caloric function defined by 

V{x, s + l) = Y e)V{x + e, s) in B2r{xQ) x {-4r^ < A; < 4r^}, 
eer 

V{x, s) = l on dp{B2r{xQ) x {-4r2 < A; < Ar'^}) n {s < 0}, 
s) = on 9p(52r(a;o) x {-4r2 < A; < 4r^}) n {s > 1}. 
By the maximum principle we get 

C/(x, s) > y(x, s), S>Q^X ^ B2r{xQ). 

Using the parabolic Harnack inequality applied to V{x^s) in B2r{xo) x 
we deduce that 

inf U{z,s) >cV{xo,-r'^) = c, 0<s<r^. 

2e-Br(a;o) 

[Note that F = 1 on B2r{xo) x {-Ar"^ <k<0}.] □ 
3.2. The global adjoint solution. 

Theorem 8. There exists a positive adjoint solution M defined globally 
in Zi'^. This solution is unique up to a multiplicative constant and verifies 

(3.10) Y M{y)<C Y ^%)' xeZ^r>0, 

where C = C{d, a, T) . 
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Proof. Let 

mi{y) = ai[Gi+i{{), y) - Gi(0, y)], y G B^, (0), / = 1,2,..., 

where G;(0, •) is the Green function of (5j)jgN in the ball i?2'(0)i ^ = 1, 2, . . . , 
with pole at the origin and where the ai are chosen so that 

(3.11) mKO) = l, / = 1,2,.... 

It is easy to see that the elhpticity condition (1.1) implies a local Harnack 
principle for the nonnegative adjoint solutions. This local Harnack principle 
and the normalization condition (3.11) imply that the mi verify 

mi{y)<C, yeB2,{0),l>k, 

with a constant C = C{k) depending only on k. The diagonal process allows 
us then to deduce the existence of a global positive adjoint solution M de- 
fined on Z"^. The fact that this global adjoint solution is unique (up to a mul- 
tiplicative constant) follows from the normalized adjoint Harnack principle 
(3.1) applied to M1/M2 where Mi and M2 denote two global positive ad- 
joint solutions. Indeed it is always possible to suppose that infgd M1/M2 = 
and associate to e > 0, G Z*^, such that Mi{ze) / M2{zs) < e. By (3.1) 
sup^^(2e) M1/M2 < Ce with a constant C > independent of R and it suf- 
fices to let i? ^ 00 and e ^ to deduce that this function is constant. Finally, 
the doubling property (3.10) follows from Theorem 7. This completes the 
proof of Theorem 8. □ 

4. The Gaussian estimates. The first step in proving the upper Gaussian 
estimate (1.5) is to prove the following mass escape estimate for (S'j)jgN 
(cf. [17, 18]). 

Lemma 2. Let {Sn)n&i be as in Section 1. Letpn{x,y), x,y £Zi'^, n = l, 
2, . . . , denote the corresponding transition kernel. Then there exist C , c> 0, 
such that 

(4.1) Pn(x,y)<Cexp -c— , x E Z"^, n, = 1, 2, . . . . 

\x-y\>R ^ " ^ 

Proof. To prove (4.1) it suffices to prove that 

(4.2) Pn{x^y)<e-'^'l^ 

l(x~y)>R 

for every linear form / : R"^ — > R such that l{x — y) < \x — y\,x,y ^ IV^. Let 
s > 0; we have 

E Pn{x,y)< E e-^^+^'(--^V„(x,y) 

l{x-y)>R l{x-y)>R 

yeZ'i 
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from which it fohows that 
Pn{x,y) 

l{x-y)>R 



(4.3) 



j/i,...,j/„-ieZ'* 
On the other hand, we have 



eer 



= l + 0{s'e^'), 

where the last equahty follows from the fact that (5'.«)„gN has symmetric 
increments. We deduce then that 



(4.4) sup 



If |s| it suffices to observe that 



(4.5) 



y'eZ'i 



<1 + Cs^ <e^' , \s\ <1. 



eer 



Putting together (4.3), (4.4) and (4.5) we deduce that 

E P„(x,y)<e-^+^-^ 

l{x-y)>R 

and optimizing over s, we deduce (4.2). □ 

The second step is to apply the parabolic adjoint Harnack principle to 



v{y,t) 



Pt{x,y) 
M{y) ' 



yGZ^t=l,2, 



This gives 



v{y,n)<C inf v{z,2n), n>C. 
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Hence (with the notation of Section 1) 

v{y,n)V{y,V^)<C ^ v{z,2n)M{z) < J2 P2n{x,z), n>C, 

in the case ^/n <c\x — y\, with c > smah enough. This imphes, in this case, 
by Lemma 2 

v(y, n)V(y, y/n) < C exp — c ) , n>C. 

\ n J 

In the case -y/n > c\x — y\ it suffices to observe that the Gaussian factor in 
(1.5) is ~ 1. Hence 

Pn{x,y) < ^^^^^ exp (-C— — —\ n>C. 



The doubhng property (3.10) allows us to symmetrize the volume factor in 
this estimate and obtain 

Let us observe that for 1 <n<C (1.5) and (1.6) are immediate consequences 
of the local Harnack estimate. This completes the proof of the upper estimate 
in Theorem 1. Finally the lower Gaussian estimate (1.6) can be deduced from 
the upper Gaussian estimate (1.5) and the parabolic adjoint Harnack by a 
standard procedure. We first use (1.5) to deduce that for ^ > large enough 

(4.6) ^ Pn{x,y)>^. 

\x—y\<Ay/n 

Parabolic adjoint Harnack applied to 

n{y,t) = P^, yez^t = i,2,..., 

M{y) 

estimate (4.6) and the doubling property (3.10) imply therefore that 

The lower off-diagonal estimate (1.6) is easily deduced from (4.7) by applying 
successively the parabolic adjoint Harnack inequality. More precisely, let us 
fix X and n as in (4.7) and let y G Z"' such that |y — x| < n/C with C > 
large enough. Let k be the smallest integer '>\x — y\'^/n. Put 

{aj,tj) = (^ttj, (^1 + J = 0, .. 
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with 

k 

Then {x^n) = {ao,n) and {y,2n) = {ak,2n). Moreover 



ao = x, ak = y, \aj+i-aj\M — - — , 0<i<A;-l. 



|aj+i - PS -, 0<j<k-l. 
Hence the parabohc adjoint Harnack inequahty yields 



n 



M{y) - M{x) - CV{x,^/^) ' \ n 
This completes the proof of Theorem 1. 

5. Proofs of results. To the process {Sj)j£]si we shall associate the cor- 
responding space-time process 

S, = {S,,to-j)eZ'^xZ, toeZ,...,i = o,l,... . 

For any cylinder Q = Q,x{a<k<b}, Q being a bounded domain in Z'^, we 
shall denote tq the first exit time of Sj from Q. The caloric measure in Q 
at (xo,to) is defined by 



(a:o,to) 

Q 



(i?)=P(,„,t„)[S,^Gi?], EcdpQ. 



Observe that for each Lp:dp{0, x {a < k < b}) — > R, the solution of the 
boundary value problem 

u{x, t + 1) = ^ 7r{x, e)u{x + e,t) in Q x {a< k <b}, 
eer 

u{x,t) = Lp{x,t) on dp{Q X {a < k <b}), 

can be represented by means of cj^O'*" = uiq°'^°\ (xo,to) € ^ x {a < k < b} 
as follows: 

n(xo,to) = ^(xo,to)['/5('5fQ)]= J2 9'(y,sK'°'*"(2/,s)- 

5.1. A lower estimate for the caloric measure. Let the notation be as in 
Section 2 and let Q = x Z. For Y = {y,s) ^ diQ, r > we shall denote 

Ar{Y) = diQnCr{Y), 

where diQ = (917 x Z. 

Lemmas. Let lo^ = Loi^'^l^.. Then 

(5.1) inf uj^(A2r(Y))>9, Y e diQ,c<r < Ro/2, 

xeQr{Y) 

where 9 = 9{d, a, F) > 0. 
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Proof. Let Y = {y,s) G diQ. It is clear that there exists a cyhnder 
C = B^r{z) X {s - 4r2 < < s} c C2r{Y) \ (0 X Z) (provided that is small 
enough). Let A' = B^j.{z) x {s — 4r^} denote the bottom of this cylinder. 
Using the maximum principle we deduce 

(5.2) ^^(A2.(y))>^;(X)=a;^^^.(^)(A') 
in Q2riy), and 

v{X)>v'{X)=oj^,{^') 
in C . On the other hand, the parabolic Harnack principle applied to v gives 
inf w^(A2r(y)) > inf viX) 

(5.3) >cv{z,s-2r'^) 

> cv {z^ s — 2r^). 

But v' can be extended from C to a large cylinder C" = B^^riz) x [s — 6r^, s] 
by 

(5.4) v'{X)=u;^„{dpC" n{t<s-Ar^}) 

so that f ' = 1 on C" n {t < s — 4r^} and the lower estimate (5.1) follows then 
from the Harnack principle. □ 

Corollary 1. Under the assumptions of Lemma 3, let u be a non- 
negative solution of Cu = in QsriY), which vanishes on A2r{Y). Then 
Mr = supQ^(y) u satisfies 

(5.5) Mr < pM2r, c<r<Ro/4, 
with a constant < p = p{d, a, F) < 1. 

Proof. Let X G QriY); we have 

n{X)= u{Z)u^^{Z)= u{Z)u''{Z). 

ZedpQ2r{Y} Z&dpQ2riY)-A2r{Y) 

Hence 

U{X) < io''[dpQ2r{Y) - A2r{Y)]M2r 
= {l-u;^'^[A2r{Y)])M2r 

< (1 - e)M2r = pM2r. □ 
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5.2. The Carleson principle for caloric functions vanishing on the bound- 
ary. Let the notation be as above. Let Y = {y,s) € diQ and c < r < Rq/2. 
Assume that u is a nonnegative caloric function in Qn (B^riv) x {s — 9r^ < 
k<s + 9r2}) and li = on diQ D {B2r{y) x {s-4r^ <k<s + Ar"^}). Then 

(5.6) u{X)<Cu{Yr), X£Qr{Y), 

where C = C{d,a,T) > 0. To prove (5.6) we first observe that the local 
Harnack principle allows us to assume that the parabolic distance of X 
from dpQ2r{Y) is sufficiently large. We shall denote by S{X) [X G C2r{Y)] 
this distance and suppose that 6{X) = Dist{X,dpQ2riY)) > C. Geometric 
considerations in combination with the parabolic Harnack principle imply 
that 

(5.7) 5^X)u{X)<Cr^u(Yr), X£Q2riY), 

where 7 and C > are positive constants depending on d, a, T. Let 
5{X) = Dist{X, dpC2riY)). Let < Eq < IQ-^ smah enough so that 

''« = (r4^<i> 

where p is the constant given by Corollary 1. We shall distinguish two cases. 
First assume that 5 = d{X) < eQ6{X). In this case we have 5 = S{X) = 
5{x,t) = dist(x, (?r2) = \x — xo\ for some xq S dQ. By Corollary 1 applied 
to u in Q2s{Xo) = C2s{xo,t) D Q, we have 

u{X) < sup u< p sup u. 

But 

6{X) < 6{Z) +46< ~5{Z) + 4eo~5{X) 
where Z = (z, r) G Q'^siXo) is such that u{Z) = supq^^(j5Cq) u. Therefore 

5{X)<{l-4e^)-^{Z) 

and this gives 

~5{Xyu{X) < (1 - 4eo)-"'p~5{Z)"'u{Z) 

(5.9) 

<9q sup 5{Xyu{X). 

X£Q2r{Y) 

It remains to examine the case where S = 6{X) > eo^(^)- In this case we 
have 

(5.10) 6{X)^u{X)<eQ^6iX)^u{X)<e^^ sup 6{Xyu{X). 

Q2r(Y) 
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Putting together (5.9) and (5.10) we deduce that 

sup 6{Xyu{X) < mayii 9q sup S'^UjEq^ sup 

Q2riY) V Q2r[Y) Q2riY) / 

Using (5.8), the fact that S{X) w r, X G Qr{y) and (5.7), we deduce the 
estimate (5.6). 

5.3. The boundary Harnack principle and proof of Theorem 4. For 
Y = {y, s) G X Z; c<r < R< Rq/2, we denote 

^rAv) = BR{y) n {x G dist(x, d^) < r}, 
DrAY) = Dn^r = ^rAv) x{s-R'<k<s}, 
^R,riY) = AR,r = dpDn^r H {x G 17, < dist(x,0r2) < r}, 

SrAY) = SR,r = dpDR^r n{x£n, dist(a;, dQ) > r}. 



Lemma 4. Let Y = {y,s) e dQ x Z. Then we have 

(5.11) Px[Srn,,^^^ e SKr,r] > Px[Srn,,„. ^ ^KrA^ X G Qr{Y), 

provided that K > Kq is large enough. 

Proof of Theorem 4. Theorem 4 is an immediate consequence of 
the Carleson principle and Lemma 4. Indeed, we may always assume that 
v{YKr) = u{Y_j^j.) = 1. By Carleson, v < cq = co{d,a,T) in QxriY) (which 
contains DktA- constant cq can be chosen so that, by Harnack, u> 1/cq 
on SKr,r- Let Mo = cqu and vq = — uq. Let X G Qr{Y). By Lemma 4 we 
have 

Vo{X) < LU^'iAKrAY)) < i^^'iSKrAy)) < M^) 

and then 

V 2 f ''^o \ 2 

sup — = Cq sup 1" 1 ^ ^Cn. 

QriY)U Qr{Y)\U0 J □ 

Proof of Lemma 4. To prove the estimate (5.11), it suffices to show 
that if n, v: DKr,r satisfy 

Cu = 0, U>0 in DKrr', U>1 on SKrr, 

(5.12) 

Cv = 0, V <1 in DKr,r\ f < on dpDKr,r \ ^Kr,r, 

then we have 

(5.13) v<u in Qr = Qr{Y) 
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provided that K > Kq is large enough. 

The first step is to prove that under (5.12), u verifies the lower estimate 

(5.14) ^(X)>25f^^^^^^l^y, X = {x,t)^Qr{Y), 



for appropriate constants 5, 7 > 0. Let y = (i?o — 5r) |^„^^| ■ We assume that 
i^>10. We define ?i:Q^^R by 

Cu = in Qer , 



u = min(ti, 1) on dpQ^r H DKr,r, 

U=l on dpQer \ DKr,r- 

Since u>0, we have, by the maximum principle, < n < 1 in Qqt, and since 
u>l on SKr,r we have u > u on Q,.- Let z G i?2r(y) satisfying 
|i- (i2o-4r)|™y| < 1. Let Z = (5, s - 4[r]2). Let w be defined by w{X) = 

l-u{X), X eU, U = QQr<^{B2r{y) X {s-4r2 < < s}). It; vanishes on diU\ 
di{B2r{y) X {s — 4r^ <k < s}). On the other hand, by the same argument as 
in Lemma 3 we see that uj^{diU \ 5/(i?2r(y) x {s — 4r^ < A; < s})) > c> 0. 
It follows then that w{Z) < OsupuW, where <0 < 1. This means that 
1 — u{Z) <6<1 and therefore u{Z) > 1 — 6 > 0. Using Harnack and the fact 
that u>u on Qr we then deduce (5.14). It follows from (5.14) that 



u{X)> 26 K-^ , XeQr\ D^^,/K , 

where we assume r > Kc and K sufficiently large. Observe that we have in 
particular u{X) > 25K~'^ , X G S^^r/K- 

The second step is to prove that there exists > such that 

(5.15) t;(X) <exp(-iVE:), X£Dr,r- 

Let J = 1, 2, . . . such that 2j + \<K and Xj such that 

Xj G dpDi2j-l)rTi sup V = v{Xj) = v{xj,Sj). 



D(2j-l)r,r 



Let U = {B2r{xj) X {sj - 8r'^ < k < sj}) n DKr,r- We have v < 
on dpU \ dp{B2r{xj) x {sj — 8r^ < A; < Sj}) and using the fact that 

u:^/ {dpU - dp{B2r{xj) x {sj - Sr"^ < k < sj})) > c> 



u 

we deduce that 



v{Xj) < sup v<9supv, 

Un{Br{xj)x{sj-4r2<k<Sj}) U 



where < ^ < 1. Hence 



sup v<6supv<p sup V, 

-D{2j-l)r,r -D{2j + l)r,r 
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where < p < 1. Iterating this estimate we obtain 

supv<p^ sup v<e~^^, 

Dr,r -D(2fe + l)r,r 

where 2k + \ < K < 2k + 'i. Thus (5.15) is proved. It fohows in particular 
that V < 6K~'^ in D^.r provided that K is large enough. 

From the previous considerations it follows that ui = > in Dj.^r/K 
and ui> 1 on Dr^r \ D^ j-jj^ (that contains Sr^r/K) = ^i'^''^ — u) ^ 

< 1 in Dr^r (that contains Dr^r/x) with < on Sr^r/K- Iii particular, 
we have 

^ N ^ 

Ul — Vl = —^{u — V) >0, on Drr\ Dj-r/K- 
' ' / 

On the other hand, ui, vi satisfy the same assumptions as u, v with r 
replaced by r/K. We can iterate and define Uj, Vj such that 

K-yy 



Uj - Vj = — j {U-V)>0 on D^/KJ^r/KJ \ Dr/K0,r/K^+^ > 

J = 1, 2, ... , and consequently 

U-V>G on 5(y) = y D^iK,^^/KJ \ D^/Kj^r/KJ+^- 

Let now Xq = {xQ,to) G Qr C Dr^r{Y) and let Xq = (xqi^o) where 
xq £ dfl satisfies dist(x,9r2) = |xo — xq\. Then DKr,r{Xo) C i^(_ft'+2)r,r(^) 
and SKr,riXo) C S(^x+2)r,r(X)- Replacing K with K + 2 in the previous con- 
siderations, we deduce that n > f on S{Xq) that contains Xq. This completes 
the proof of Lemma 4. □ 

5.4. The boundary backward Harnack principle. Let the notation be the 
same as in Theorem 5. Let Br{yo) = Br- First we observe that the Carleson 
principle in combination with the parabolic Harnack principle give 



(5.16) ui X 



r^2 



< C min n, x £ Br — B^/^- 

fir/2 x{''-V4<t<8r2} ' 



On the other hand, by Harnack 

(5.17) max u<C min u. 

Br/2X{i=[r2/8]} Br/2x{r2/4<t<8r2} 

Since u = on dBr x N, by (5.16), (5.17) and the maximum principle we 

get 

max u<C min u. 

BrX{ [r2 /8] <t<8r2} 3^/2 x {■r2/4<t<8r2 } 



xo £ (x, t) £ Br X {r"^ < t < Sr^}. 
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In particular, we have 

(5.18) max u<C min u. 

■Br/2X{[rV8]<t<8r2} B,./2 X {r2/4<t<8r2} 

Let now v be another nonnegative caloric function in x N such that v 
vanishes on the lateral boundary dBr x N. Then, there exists C > such 
that 

u{xo, [rf)v{x,t) < Cv{xoA[rf)u{x,t), 

(5.19) 

"r/2 

To prove (5.19) we first use a covering argument to cover di{Br x {r^ <t< 
3r^}) with cylinders Cj = B^riUj) x {sj — e^r^ <t < Sj}, j = 1, . . . , N , where 
e > is chosen sufficiently small and where Y^^^ = {yj,Sj) € di{Br x {r^ < 
t < 3r^}), and apply in each of these Cj the boundary Harnack principle to 
get 

(5.20) uixl^f^Mx, t) < Cv{y['^,)u{x, t), (x, t) G Cj, 

where Yj- , Y^ , r > 0, are defined as in Section 2. By Harnack, we have 

(5.21) v{Y^^/,)<Cv{xoA[r?), 

(5.22) u{Yl%)>cu{xo,[r^/A]). 
Using (5.18) we deduce from (5.22) that 

(5.23) u{Yl%)>cu{xo,[r?)- 
Putting together (5.20), (5.21) and (5.23) we deduce that 

u{xo, [rf)v{x,t) < Cv{xo,A[rf)u{x,t), 

(5.24) 

Xo G Br/2, {x,t) £ Cj,j = 1,...,N. 

On the other hand, we have, by Harnack, 

(5.25) v{x,t) <Cv{xo,i[r]'^), r'^<t<3r'^, dist{x, dBr) > 6r, 
where < 6 < 1/2. Again, by Harnack combined with (5.18) 

(5.26) u{x,t)>cu{xo,[r]'^), r'^<t<3r'^, dist{x, dBr) > Sr 

and (5.19) follows from (5.24)-(5.26) with an appropriate choice of (5 > 0. 
We are now able to get the estimate (2.3). We shall use (5.18), (5.19) and a 
time-shifting argument. Let u be as in (2.3) and let v{x,t) = u{x,t + 2[r]^). 
By (5.19) we have 

u{xo, [rf)u{x, t + 2[rf) < Cu{xo, d[r]'^)u{x, t), 

(5.27) 

{x,t) £ Br X {r2 < t < 3r2}, 
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with xq G Bj,/2 fixed. Equation (2.3) follows from (5.27) and the estimate 

u{xo,6[rf)<Cu{xo,[r?), 

which is an immediate consequence of (5.18). This completes the proof of 
Theorem 5. 
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